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Abstract. We outline the key-steps toward the construction of a physical, fully 
relativistic cosmology. The influence of inhomogeneities on the effective evolution 
history of the Universe is encoded in backreaction terms and expressed through 
spatially averaged geometrical invariants. These are absent and potential candidates 
for the missing dark sources in the standard model. Since they can be interpreted as 
energies of an emerging scalar field (the morphon) , we arc in the position to propose a 
strategy of how phenomcnological scalar field models for Dark Energy, Dark Matter and 
Inflation, that are usually added as fundamental sources to a homogeneous-geometry 
(FLRW) cosmology, can be potentially traced back to inhomogeneous geometrical 
properties of space and its embedding into spacetime. We lay down a line of arguments 
that is - thus far only qualitatively - conclusive, and we address open problems of 
quantitative nature, related to the interpretation of observations. 

We discuss within a covariant framework (i) the foliation problem and invariant 
definitions of backreaction effects; (ii) the background problem and the notion 
of an effective cosmology; (iii) generalizations of the cosmological principle and 
generalizations of the cosmological equations; (iv) dark energies as energies of an 
effective scalar field; (v) the global gravitational instability of the standard model 
and basins of attraction for effective states; (vi) multiscale cosmological models 
and volume acceleration; (vii) effective metrics and strategies for effective distance 
measurements on the light cone, including observational predictions; (viii) examples of 
non-perturbative models including explicit backreaction models for the LTB solution, 
extrapolations of the relativistic Lagrangian perturbation theory, and scalar metric 
inhomogeneities. The role of scalar metric perturbations is critically examined and 
embedded into the non-perturbative framework. 
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1. General relativity and cosmology 

1.1. The foliation issue and the notion of an effective cosmology 

The homogeneous-isotropic standard model of cosmology, being itself a particular 
solution of Einstein's general theory of relativity, does by far not exploit the degrees 
of freedom inherent in the geometry as a dynamical variable. It is this richer tone of 
general relativity - as compared to the Newtonian theory - that opens the possibility 
to generalize cosmological models, notably by including inhomogeneous structure also 
in the geometrical variables. There are several guidelines to be emphasized in such 
a generalization: firstly, a cosmology is thought of as an evolving space section that 
implies the need to speak of a foliated spacetime, introducing four degrees of freedom 
(the lapse and shift functions in an ADM setting). This dependence on the foliation 
should not be confused with the coordinate- or gauge-dependence of the resulting 
cosmological equations and variables, however. Secondly, a cosmology purports an 
effective point of view in the sense that the evolving spatially inhomogeneous variables 
are thought of as being "averaged over" in a way that has to be specified. We aim 
at a description that only implicitly refers to a metric. However, if a metric is to be 
specified, a cosmological metric is then to be considered as an effective, "smoothed 
out" or template metric, being not necessarily a solution of the equations of general 
relativity. Finally, a physical cosmology should be characterized by such an effective 
evolution model, an effective metric to provide the distance scale for the interpretation 
of observations, or alternatively an evolution model for average characteristics on the 
light cone, together with a set of constraint initial data. These latter are to be related 
to physical properties of fundamental sources, but also to the geometrical data at some 
initial time (effective, i.e "averaged" quantities of known energy sources, but also of 
intrinsic and extrinsic curvature). This latter clearly emphasizes the absence of any 
phenomenological parameters. Those would just parametrize our physical ignorance. 
All these points will be made explicit in what follows. 

1.2. The dark side of the standard model: postulated sources and proposed solutions 

The high level of idealization of the geometrical properties of space in the standard 
model leads to the need of postulating sources that would generate "on average" a 
strictly, i.e. globally and locally, homogeneous geometry. It is here where a considerable 
price has to be paid for a model geometry that obviously is not enough to meet physical 
reality. Assuming a FLRW geometry 96 percent of the energy content is missing in 
the form of a) a postulated source acting attractive like matter, so-called Dark Matter 
(= 23 percent) and b) a postulated source acting repulsive, so-called Dark Energy 
(= 73 percent). Evidence for the former does indeed come from various scales (galaxy 
halos, clusters and cosmological, see e.g. [Ill] ), while evidence for the latter only comes 
from the apparent magnitude of distant supernovae (see [HU EHl [60] for the latest data) 
that, if interpreted within standard model distances, would need an accelerating model. 
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In the simplest case this volume acceration is achieved by a homogeneous-isotropic 
cosmology with a cosmological constant. It should be emphasized that when we speak 
of evidence, we already approach this evidence with model priors [71 1 11161 1123] . Keeping 
this idealization for the geometry of the cosmological model on all scales, one has to 
conjecture fundamental fields e.g. in proportion to the missing Dark Matter within 
small-scale systems also on cosmological scales. The search for these fields is one major 
research direction in modern cosmology. 

Another huge effort is directed toward a generalization of the underlying theory 
of gravitation. While this would generalize the geometry of the model, it is not clear 
why most of these efforts go into a generalization of general relativity and not into the 
generalization of the cosmological model within general relativity. There are certainly 
good lines of arguments and various motivations in particle physics and quantum gravity 
to go beyond the theory of Einstein (for reviews see [19], |114] ). but the "dark problem" 
may be first a classical one. 

Looking at generalizations of the standard model within general relativity can be 
identified as a third research direction to which we dedicate our attention here. In 
light of current efforts it is to be considered conservative, since it does not postulate 
new fundamental fields and it does not abandon a well-tested theory of gravitation 
[28] , [101] , [77] (for reviews on the physical basis of this third approach see [TTJ [23] and 
[102] ). Among the works in this latter field, research that analyzes spherically symmetric 
exact solutions has been meanwhile developed to some depth, and has determined 
the constraints that are necessary to explain Dark Energy, on a postulated observer's 
position within a large-scale void (see [581 ESI El I117[ 1971 193] and references therein, 
as well as the contributions [9] and [89] in this volume). 

1.3. Fictitious and physical backgrounds: a more realistic cosmological principle 

Perhaps a reason for not questioning the standard model geometry within general 
relativity and to go for the search for fundamental fields or for generalizations of the 
laws of gravitation is the following belief: effectively, i.e. "on average", the model 
geometry has to be homogeneous, since structures should be "averaged over". Then, 
due to observational facts on large scales (the high degree of isotropy of the Cosmic 
Microwave Background, if the dipole is completely eliminated due to our proper motion 
with respect to an idealized exactly isotropic light sphere), and first principle priors (the 
strong cosmological principle that requires the universe model to look the same at all 
places and in all directions), the model's geometry - supposed to describe the Universe 
on average - is taken to be locally isotropic. 

Taking this reasoning at face value we must note the following. The notions of 
homogeneity and isotropy in the standard model are both too strong to be realistic: 
firstly, local isotropy implies a model that is locally and globally homogeneous, i.e. 
despite the conjecture that the homogeneous model describes the inhomogeneous 
Universe "on average" , this strict homogeneity does not account for the fact that any 
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averaging procedure, in one way or another, would introduce a scale-dependence of the 
averaged (homogeneous) variables [57]. This scale-dependence, inherent in any physical 
averages, is suppressed. Even if a large scale of homogeneity exists (we may call this 
weak homogeneity principle), the model is in general scale-dependent on scales below 
this homogeneity scale |120j (see the contribution by Sylos Labini in this volume |121j ). 
The same is true for isotropy: while the averaged model may be highly isotropic on 
large scales, a realistic average distribution on smaller scales is certainly not (we may 
call this weak isotropy principle). Correspondingly, a weak cosmological principle that 
applies only to the largest scales would be enough to cover the reality needs while still 
respecting observational evidence. 

We may summarize the above thoughts by noting that, on large scales, a 
homogeneous-(almost)isotropic state does not necessarily correspond to a homogeneous- 
(almost)isotropic solution of Einstein's equations. These former states are the averages 
over fluctuating fields and it is only to be expected that the state coincides with a strictly 
homogeneous solution in the case of absence of fluctuations. In other words, looking at 
fluctuations first requires to establish the average distribution. Only then the notion 
of a background makes physical sense [THl [79] ; see also the contribution by Kolb [80] in 
this volume. Current cosmological structure formation models, perturbation theories or 
N-body simulations, are constructed such that the average vanishes on the background 
of a homogeneous-isotropic solution [I6]. A such chosen reference background may be 
a fictitious background, since it arises by construction rather than derivation. On the 
contrary, a physical background is one that corresponds to the average (whose technical 
implementation has to be specified, and which is nontrivial if tensorial quantities like the 
geometry have to be "averaged" ). A sound implementation of a physical background will 
be a statistical background where not only solutions but also the ensembles of solutions 
are averaged. Having specified such an averaging procedure, a physical cosmological 
model may then be defined as an evolution model for the average distribution. Despite 
these remarks it is of course possible that the homogeneous solution forms at the same 
time the average. A well-known example is Newtonian cosmology [I6]. It is also 
conceivable that a Friedmannian background provides, in some spatial and temporal 
regimes, a good approximation for the average (compare here the analysis of the stability 
properties of Friedmannian backgrounds |113j in this volume). 

2. Refurnishing the cosmological framework 

2.1. Effective evolution of inhomogeneous universe models 

Taking the point of view of generalizing the cosmological model within general relativity 
by abandoning the strong cosmological principle (strict homogeneity and isotropy 
on all scales) and replacing it by the weak cosmological principle (existence of a 
homogeneity scale and restriction to effective states that are almost isotropic on the 
scale of homogeneity) leads us to a "rewriting of the rules" to build the cosmological 
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model. We shall consider the rules that led to the standard model of cosmology 
and replace them by their more general counterparts. It follows a basically similar 
framework that displays, however, a signature of inhomogeneity through the occurence 
of so-called backreaction terms and through a manifest scale-dependence. We shall not 
introduce new principles or assumptions, apart from the above outlined relaxation of the 
cosmological principle. We shall restrict ourselves to the simplest case of an irrotational 
dust model, except in the last section of this review (for generalizations of the dust model 
[18] with non-constant lapse function see [I9], and for additionally non-vanishing shift 
see [T0l[Tn[83lH2l[T22l[63l[88]l 

2.2. The key-steps for generalizing the cosmological framework 

We shall now paraphrase the foundations of the cosmological equations of the standard 
model and give their generalized counterparts. 

• As in the standard model we introduce a foliation of spacetime into fiow-orthogonal 
hyp ersurf aces. We generalize the notion of Fundamental Observers to those that are in 
free fall also in the general spacetime. Although, as in the standard model, this setup 
depends on the chosen foliation, we presume that this choice is unique as it prefers the 
fundamental observers against observers that may be accelerated with respect to the 
hypersurfaces. A general inhomogeneous hypersurface - contrary to the homogeneous 
case - will, in this setting, unavoidably run into singularities in the course of evolution. 
This is to be expected in a given range of spatial and temporal scales, since we are 
treating the matter model as dust. This is not a problem of the chosen foliation, but 
a problem of the matter model that has to be generalized, if small-scale structure 
formation has to remain regular, and this can be achieved by the inclusion of velocity 
dispersion and vorticity. 

• As in the standard model we confine ourselves to scalar quantities. We replace, 
however, the homogeneous quantities by their spatial averages, e.g. the homogeneous 
density Qnit) is replaced by (f?)^ (t) for the inhomogeneous density g that is volume- 
averaged over some compact domain V. We realize the averaging operation by a 
mass-preserving Riemannian volume average. In some mathematical disciplines and 
in statistical averages at one instant of time, it may be more convenient to introduce 
a volume-preserving averager, but thinking of an averaging domain that is as large 
as the homogeneity scale we have to preserve mass rather than volume. Furthermore, 
the average is performed with respect to the above-defined Fundamental Observers. 
Spatially averaging a scalar X*), as a function of Gaussian coordinates and a 
synchronizing time t, is defined as: 

(vl/(t,X'))^(t):=-l / ^{t,X^)dfi, , (1) 

with the Riemannian volume element dfig := ^/gd^X, g := det{gij), and the volume of 
an arbitrary compact domain, := j^^/gd^X. Note that within a more general 

setup that includes lapse and shift functions, we would have to consider the question 
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whether the locally lapsed time is replaced by a global "averaged time" that would 
involve an average over the lapse function |109j . Here, the dust cosmology is already 
synchronous, so that this question does not arise. Note furthermore, that the building 
of averages is done in the inhomogeneous geometry. The averages functionally depend 
on the inhomogeneous metric, but this latter needs not to be specified. We may talk of 
a kinematical averaging that does not deform the geometry, i.e. that does not change 
the physical properties of the inhomogeneous spacetime. For other strategies, see |57] , 
and references therein. Section |6l as well as the contribution of Wiltshire |130] in this 
volume. 

• We generalize the kinematical laws of the standard model a) for the volume expansion 
(the Hamiltonian constraint in the ADM formulation of general relativity) and b) 
for the volume acceleration (Raychaudhuri's equation in the ADM formulation of 
general relativity) by dropping the symmetry assumption of local isotropy. The 
general equations are then volume-averaged, leading to the following general volume 
expansion and volume acceleration laws (for a volume scale factor, defined by a-ri (t) := 
iyvit)/Vx>{ti))^^^; the overdot denotes partial time-derivative, which is the covariant 
time-derivative here) [T8] : 

3^ = -4nG{g)^ + Qv + A ; 3H^ + ^ = 8nG {g)^ - ^Wv ~ ^Qv + A , (2) 

where Hj) denotes the domain dependent Hubble rate i^x" = cix>/ax> = —1/3 {K)jy, 
K is the trace of the extrinsic curvature Kij of the embedding of the hypersurfaces 
into the spacetime, and A the cosmological constant. The kinematical hackreaction 
Qt> is composed of averaged extrinsic curvature invariants, while W-p is an averaged 
intrinsic curvature invariant that describes the deviation of the average of the full (three- 
dimensional) Ricci scalar curvature TZ from a constant-curvature model, 

Qv:={K'-K)K\)^-l{K)l ; := (7^)^ - ^ . (3) 

The kinematical backreaction Qx> can also be expressed in terms of kinematical 
invariants, where the extrinsic curvature is interpreted actively in terms of (minus) 
the expansion tensor: 

Qv:=lmv-{e)l)-2{aX, (4) 

where 6 is the local expansion rate and cr^ := l/2aija^^ is the squared rate of shear. Note 
that Hx) is now defined as if-p = 1/3 (^)p,. Qv appears as a competition term between 
the averaged variance of the local expansion rates, {6'^)x>~ {^)'t>^ the averaged square 
of the shear scalar ((T^)^ on the domain under consideration. 

For a homogeneous domain the above backreaction terms Qt> and Wd, being 
covariantly defined and gauge invariants (to second order) in a perturbation theory 
on a homogeneous background solution, are zero. They encode the departure from 
homogeneity in a coordinate-independent way [851 IMl EH] • 
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The integrability conditions connecting the two Eqs. ([2]), assuring that the 
expansion law is the integral of the acceleration law, read: 

{g)v + ^Hv{g)^ =0 ; a^\alWv) + a^^a^Qv) = 0. (5) 

While the mass conservation law for the dust is sufficient in the homogeneous case, there 
is a further equation connecting averaged intrinsic and extrinsic curvature invariants 
in the inhomogeneous case. The expressions in brackets are conformal invariants (for 
further details see [23]). 

2.3. Interpretation and Discussion 

The interpretation of these average equations as generalized or evolving backgrounds 
[23] , [79] implies that the new second conservation law in Eq. ([5]) describes an 
interaction between structure formation and background curvature. In the standard 
model this latter is absent and structures evolve independently of the background 
having homogeneous geometry. This homogeneous curvature background furnishes 
the only solution of ([5]), in which structure formation decouples from the background 
(the expressions in brackets in the second conservation law are separately constant). 
Backreaction on such a fixed background decays in proportion to the square of the 
density and is unimportant in the Late Universe [HI [201 123]. This degenerate case of a 
decoupled evolution explains the fact that in Newtonian and quasi-Newtonian models 
backreaction has no or little relevance [23]; in the Newtonian case [16], as well as in 
Newtonian [251 176] and spatially fiat, relativistic spherically symmetric dust solutions 
(see Section [7]), Qx> vanishes. In models with homogeneous geometry and with periodic 
boundary conditions imposed on the inhomogeneities on some scale, the backreaction 
term is globally zero and describes cosmic variance of the kinematical properties in the 
interior of the periodic universe model. 

In general, a physical background "talks" with the fiuctuations, and it is this 
coupling that gives rise to an instability of the constant-curvature backgrounds as we 
discuss below. The essential effect of backreaction models is not a large magnitude of 
Qti, but a dynamical coupling of a nonvanishing Qd to the averaged scalar curvature 
deviation W©- This implies that the temporal behavior of the averaged curvature 
deviates from the behavior of a constant-curvature model. In concrete studies, as 
discussed in Section |5l this turns out to be the major effect of backreaction, since 
it does not only change the kinematical properties of the cosmological model, but also 
the interpretation of observational data as we explain in Section [61 

3. Scalar field language for backreaction: the morphon 

3.1. Rewriting the averaged equations as an effective Friedmannian model 

We rewrite the above set of spatially averaged equations together with their integrability 
conditions by appealing to the kinematical equations of the standard model, which will 
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now be sourced by an effective perfect fluid energy-momentum tensor 



(6) 



where the effective densities are deflned as 

gf,:={g)^ + g^ ; - -^^2^ " J^^^ 5 

In this form the effective equations suggest themselves to interpret the extra fluctuating 
sources in terms of a scalar fleld [T71 [211 EI]) which refers to the inhomogeneities in 
geometrical variables. 

3.2. Scalar field emerging from geometrical inhomogeneities 

In making this suggested analogy concrete we, thus, choose to consider the averaged 
model as a (scale-dependent) "standard model" with matter source evolving in a mean 
field of backreaction sources. This mean fleld we call the morphon field, since it captures 
the morphological (integral-geometrical [23]) signature of structure. (Note that in more 
general cases, involving lapse and shift functions, the structure of the scalar fleld theory 
suggested by the equations may no longer be a minimally coupled one.) We rewrite |31j : 

= + Uv ; P^ = eU>l - , (8) 

where e = +1 for a standard scalar fleld (with positive kinetic energy), and e = — 1 
for a phantom scalar fleld (with negative kinetic energy; if e is negative, a "ghost" can 
formally arise on the level of an effective scalar fleld, although the underlying theory 
does not contain one; note also that there is no violation of energy conditions, since 
we have only dust matter). Thus, from the above equations, we obtain the following 
correspondence that can be employed to change between the languages: 

-^Qv = e^l-Uv ; --i-Wp = 3f/p . (9) 

Correspondence ([9]) recasts the integrability conditions ([5]) into a (scale-dependent) 
Klein-Gordon equation for $x>; and <I>x) 7^ 0: 

$2, + 3i/^$p + e-|-t/($i,, (f?)p) = . (10) 

We appreciate that the deviation of the averaged scalar curvature from a constant- 
curvature model is directly proportional to the potential energy density of the 
scalar fleld. Averaged universe models obeying this set of equations follow, thus, a 
Priedmannian kinematics with a fundamental matter source, and an effective scalar fleld 
source that reflects the shape of spatial hypersurfaces and the shape of their embedding 
into spacetime. Given the potential in terms of the variables of the averaged system, 
the evolution of these models is flxed (the governing equations are closed). This also 
potentially flxes coupling parameters, since all involved flelds can be traced back to the 
initial value problem of general relativity. 
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3.3. Interpretation and Discussion 

The morphon formulation of the backreaction problem offers a good interpretation in 
terms of energies: a homogeneous model, Qv = (a necessary and sufficient condition 
to also drop the scale-dependence, if required on every scale), is characterized by the 
virial equilibrium condition: 

2EZ^ + El, = --^ , = ; EZ^ = e^^Vv , E^, = -UvVj, . (11) 

Deviations from homogeneity, Qx> 7^ 0, thus invoke a non-equilibrium dynamics of 
the morphon in its potential that is dictated by the effective intrinsic curvature of the 
space in which the fluctuations evolve. The off-equilibrium state can be measured by a 
Kullback-Leibler distance [TQl [2] , an entropy functional that arises naturally from the 
non-commutativity of averaging and the density evolution [7D]. It is conjectured [7D] 
and shown to hold in some popular models [95] that this entropy grows in the course of 
structure formation. 

Morphon energies are redistributed and can be assigned to the dark energies. 
Dependent on the signs of the backreaction terms (and a sign change may occur 
in the course of structure formation and by looking at different spatial scales) the 
morphon can act as a scalar field model for Dark Matter, a quintessence model for 
Dark Energy, or it may even play the role of a classical inflaton |33]. (For the different 
interpretations of scalar fields see the review |19], and for unified views the selection of 
papers ^ HOOl UlSl [m], and for Scalar Dark Matter e.g. [9H H [92]). 



4. Global gravitational instability of the standard model backgrounds 

4.1. The phase space of exact background states 

The space of possible states of an averaged cosmological model, or the space of 
physical backgrounds has one dimension more than the space of possible homogeneous- 
isotropic solutions in the standard model framework. This can be seen by introducing 
adimensional "cosmological parameters". We divide the volume-averaged expansion 
law by the squared volume Hubble functional Hx> := ax>/at> introduced before. Then, 
the expansion law can be expressed as a sum of adimensional average characteristics: 

nZ + + + + = , with : (12) 

P _ 87rG(g)p . ._ A . ^2, fep, . ._ _ . ._ _ (^',^ 

Taking for simplicity the constant-curvature parameter and the curvature deviation 
into a single full curvature parameter, + =: the generalized model offers a 
cosmic quartet of parameters. Furthermore, if we put A = 0, the expansion law defines, 
for each scale, a two-dimensional phase space of states. A one-dimensional subset of this 
phase space is formed by backgrounds with Friedmannian dynamics (for illustrations see 
Figure 3 in [31] or Figures 1 and 2 in [23], and especially Figures 1 and 2 in |113] ). 
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We can analyze the fix points and their stabihty properties in the general dynamical 
system [31]; see the detailed investigation in this volume |113] . Corresponding results for 
the LTB solutions support these findings [118] and |119j . The principal outcome of these 
studies is that the standard zero-curvature model forms a saddle point; of particular 
interest are two instability sectors for the standard model, regarded as averaged state: 
firstly, perturbed homogeneous states are driven into a sector of highly isotropic, 
negative curvature and accelerated expanding "backgrounds" where backreaction thus 
mimics Dark Energy behavior over the domain V; secondly, perturbed homogeneous 
states are driven into a sector of highly anisotropic, positive curvature, collapsing and 
decelerated "backgrounds" where backreaction thus mimics Dark Matter behavior over 
the domain V. Concrete models, discussed in Section |5] show that the former happens 
on large scales, and the latter on the scales of galaxy surveys, and also on smaller scales; 
(compare here with the phase space orbits depicted in Figure 7 in |125] ). 

Thus, qualitatively, the instability sectors identified comply with the aim to trace 
the dark components back to geometrical properties of space, but they also agree 
with the expected properties of the structure: isotropic, accelerating states on large 
scales, and highly anisotropic structures on the filamentary distribution of superclusters. 
Moreover, the curvature properties also meet the expectations (to be inferred from the 
averaged Hamiltonian constraint): on large scales the Universe is void-dominated and, 
hence, dominated by negative curvature, while on intermediate scales over-densities are 
more abundant and are individually characterized by positive curvature. (A more refined 
classification of instability sectors, associated with the dark sectors of the concordance 
model, is provided in [113] in this volume.) 

4-2. Dark Energy and Dark Matter hidden in the geometry of space 

The fact that the standard model can be globally unstable in the phase space of averaged 
states, and the fact that the instability sectors lie in the right corners to explain Dark 
Energy and Dark Matter behavior, are both strong qualitative arguments to expect that 
the conservative explanation of the dark energies through morphon energies is valuable. 
The underlying mechanism is indeed based on the fundamental existence of the relation 
between geometrical curvatures and sources dictated by Einstein's equations. 

Whether this mechanism is sufficient in a quantitative sense is to date still an open 
issue. The difficulty to construct quantitative models is to be seen in the need for 
non-standard tools, for example perturbation theory on a fixed reference background 
should be replaced by a fluctuation theory on an evolving background that captures the 
average over the fluctuations. The question whether perturbations are small can only 
be answered if we know with respect to which background they are small. Furthermore, 
since backreaction affects the geometry, it will change the interpretation of observational 
data, a problem that is intimately related to the generalization of the cosmological 
model, and to which we shall come in Section [6l 

Before, we shall in the next section explain the identifled mechanism by discussing 
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some physical properties of structure formation and its relation to the interpretation of 
geometrical curvature invariants, and how exactly they mimic the dark sources. We here 
touch on a deeper problem: backreaction effects qualitatively mimic both, Dark Energy 
and Dark Matter, simultaneously. Whether, on a given domain, or on an ensemble of 
domains on a given scale, the morphon mimics Dark Energy or Dark Matter behavior, 
changes as a function of time and as a function of scale. Moreover, the small-scale 
contribution to e.g. a Dark Matter behavior requires more sophisticated relativistic 
models than the dust model used throughout here (e.g. [Ml IllOl 138]). Considering 
rotation curves of galaxy halos, missing gravitational sources in clusters or missing 
sources on cosmological scales always needs different modeling strategies. We try in the 
following to provide a first step of disentangling Dark Energy and Dark Matter behavior 
by explicitly constructing an effective multiscale cosmological model. 

5. Multiscale cosmology: generic volume partitioning of the Universe 

5.1. A note on the non-local nature of averaging 

Contrary to the standard model, where a homogeneous background is used as 
a standard of reference for the expansion history of the Universe, a background 
constructed as the average over fluctuating fields introduces a subtle element: while 
a homogeneous geometry can be characterized locally, an average is non-local, since 
it is determined by the inhomogeneities inside, but also outside the averaging domain, 
reflecting the non-local nature of gravitation. Furthermore, an average incorporates 
correlations/fluctuations of the local fields, expressed e.g. through the variance between 
local and averaged quantities. It is this latter which is the key-driver of a repulsive 
effective pressure that arises in the averaged models, as we explain now. 

5.2. Structure-emerging volume acceleration 

The simple fact that the local expansion rates differ from their average value on 
some scale provides the reason why backreaction can produce a volume-accelerating 
component despite the decelerating nature of the general local acceleration law. We 
here are not talking about an exotic ingredient that produces such a repulsion; it is a 
basic physical property of a lumpy matter distribution as was already noted in [17] . 

This physical property can easily be understood by comparing the local and the 
volume-averaged Raychaudhuri equation (for vanishing vorticity and pressure that both 
would also act accelerating on the local level, but only on small scales) [18] : 

^ = A - AnGg + 211 -f ; (^) = A - AnG {q)^ + 2 (11)^ - (I) J , (14) 

where we defined the principal scalar invariants of the expansion tensor 6jj, 211 : = 
2/3^^ — 2cr^ and I := 9. Clearly, by shrinking the domain to a point, both equations 
agree. However, evaluating the local and averaged invariants. 
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2 (11)^ - (I) 



7? ' 



(15) 



gives rise to two additional, positive-definite fluctuation terms, where that for the 
averaged expansion variance enters with a positive sign. Thus, the time-derivative 
of a (on some spatial domain V) averaged expansion may be positive despite the fact 
that the time-derivative of the expansion at all points in V is negative. 

In concrete models this variance is the source of a possible large-scale volume- 
acceleration that would be assigned to Dark Energy in the standard model, while the 
averaged shear variance mimics an attractive source that would be missing as Dark 
Matter in the standard model on cosmological scales. Both terms are competing in the 
backreaction term Qd. Since the latter depends on scale, it may act in both ways. 

5.3. Scale-dependence made explicit 

We can go one step further and make the scale-dependence explicit by introducing a 
union of disjoint over-dense regions Ai and a union of disjoint under-dense regions S, 
which both make up the total (homogeneity-scale) region V. The averaged equations 
can be split accordingly yielding for the kinematical backreaction [29] : 



where A;^ := \ / \V\ denotes the volume-fraction of the over-dense regions compared 
with the volume of the region T>. In a Gaussian random fleld this fraction would be 0.5 
and would gradually drop in a typical structure formation scenario that clumps matter 
into small volumes and that features voids that gradually dominate the volume in the 
course of structure formation. 

Ignoring for simplicity the individual backreaction terms on the partitioned 
domains, the total backreaction features a positive-deflnite term that describes the 
variance between the different expansion histories of over- and under-dense regions. 
It is this term that generates a Dark Energy behavior over the domain V |125] (see 
also for a model by Rasanen, and p^H [T28l [T29l EM Wiltshire's model that 
is based on this term only, but includes a phenomenlogical lapse function to account 
for different histories in Ai and S regions that, this latter, we cannot implement in the 
synchronous foliation of a multiscale dust model; the reader may find more details in 
the contributions by Rasanen |108] and Wiltshire |130] in this volume). If we model 
non-zero individual backreaction terms by an extrapolation of the leading perturbative 
mode in second-order perturbation theory [851 ESj , which also corresponds to the leading 
order in a Newtonian non-perturbative model [25] , then we even produce a cosmological 
constant behavior over the homogeneity scale V, see Figure 3 in [125] . In other words, 
the fact that, physically, over-dense regions tend to be gravitationally bound, i.e. do 
not partake significantly in the global expansion, together with a volume-dominance 
of under-dense regions, already produces a large-scale kinematical pressure as a source 
of volume acceleration. A homogeneous background simply cannot account for this 
physical property. 



Qv 



^mQm + (1 - A^) Qe + 6A^ (1 - A^) [Hm - Hef , 



(16) 
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6. Effective metrics and light cone distances 

6.1. Template metrics and effective distances 

For the construction of an effective cosmological evolution model, as outlined above, a 
metric needs not be specified. The need for the construction of an effective metric in 
these models arises, since measured redshifts have to be interpreted in terms of distances 
along the light cone. Given an explicit, generic and realistic, inhomogeneous metric, the 
need for the construction of effective metrics does not arise. Also, if we succeed to 
understand the evolution of light cone averages in relation to distances, then also here 
an explicit metric will not be needed |105[I106] . Work is in progress to construct effective 
equations on light fronts by surface averaging optical scalars [36]. In this line Gasperini 
et al. in a recent paper made a first step toward defining suitable covariant and gauge 
invariant light cone averages |64] . 

The idea of an effective cosmological metric comes from the "fitting problem", 
that has been particularly emphasized by George Ellis already in the 70's [5l]. The 
observation was that an inhomogeneous metric does not average out to a homogeneous 
metric that forms a solution of general relativity. Not only the nonlinearity of the theory, 
but also simple arguments of a non-commutativity [57] between evolution equations 
and the averaging operation, give rise to the need to find a "best-fit", we may call 
it "template" geometry, that inherits homogeneity and (almost-) isotropy on the large 
scales and, at the same time, incorporates the inhomogeneous structure "on average" 
(for earlier practical implementations of this problem see [56], [M], [61], [71], [37], [62] ; 
compare here the introduction to early work on the backreaction problem, to the fitting 
problem and the discussion on geometrical optics by George Ellis in this volume |55j). 

For the solution of the fitting problem various strategies have been proposed (see [57] 
and references therein). One strategy, that allows to explicitly perform a "smoothing" 
of an inhomogeneous metric into a constant-curvature metric at one instant of time, 
is based on Ricci-fiow theory: one notices that a smoothing operation of metrical 
properties can be put into practice by a reseating of the metric in the direction of 
its Ricci curvature. The scaling equations for realizing this are well-studied, and the 
rescaling flow results in a constant-curvature metric that carries "dressed" cosmological 
variables [27], [28]. These incorporate intrinsic curvature backreaction terms describing 
the difference to the "bare" cosmological parameters as they are obtained through 
kinematical averaging. 

6.2. Reinterpretation of observational data 

The standard method of interpreting observations is to construct the light cone ds"^ = 
from the line-element ds"^ = —dt"^ + gf°'^dX^dX\ where the coefficients gf""^ are given in 
the form of a constant-curvature (FLRW) metric, and then to calculate the luminosity 
distance dL{z) in this metric for a given observed redshift z. Assuming this metric 
for the inhomogeneous Universe implies the conjecture that the FLRW metric is the 
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correct "template" of an effective cosmological metric. However, the integrated exact 
equations (the integral properties of a general inhomogeneous model) are not compatible 
with this metric, simply because the averaged curvature is assumed to be of the form 
(7^)p = G/ca"^ on all scales. Improving the metric template slightly, by replacing 
the global scale factor a{t) through the volume-scale factor aT>{t) and the integration 
constant k through the domain-dependent integration constant /cx>., one renders this 
metric implicitly scale-dependent [99]. As we explained, this is not enough since the 
averaged curvature couples to the inhomogeneities and in general deviates from the a^^- 
behavior. What we can do as a first approximation, and this would render the metric 
compatible with the kinematical average properties, is to introduce the exact averaged 
curvature in place of the constant curvature in this metric form [82]. 

The resulting effective spacetime metric consists of a synchronous foliation of 
constant-curvature metrics that are, however, parametrized by the exact integral 
properties of the inhomogeneous curvature, thus they "repair" the standard template 
metric as for the evolution properties of spatial variables. Such a construction can 
be motivated by Ricci-flow smoothing, that guarantees the existence of smoothed-out 
constant curvature sections at one instant of time, and by assuming that the intrinsic 
backreaction terms are subdominant, so that we can parametrize the metric by "bare" 
kinematical averages. To stack these hypersurfaces together introduces, however, an 
inhomogeneous light cone structure [9l], |103] . Ideally, one would wish to smooth the 
light cone too, which is also possible by employing Ricci flow techniques. Improving 
this first approach to a template metric is needed and this is work in progress. 

The result of employing improved template metrics as described above is a change 
in the luminosity distance. It will, e.g., take care of the fact that light mostly propagates 
in under-dense regions of negative curvature. This will alter the interpretation of all 
observational data formerly based on FLRW distances. A comparison of the luminosity 
distances in the multiscale models investigated in |125] . that are based on the template 
metric of [82] , with a flat ACDM model is presented in Figure 2 of |125j , see also |126] . 
Clarkson's C-function [H] (see the review by George Ellis in this volume [55]) features 
a clear minimum at redshifts of around 3 — 5, which may serve as observational evidence 
for the effective cosmologies, as proposed in [82] . As the investigated multiscale models 
show, it is not even necessary to measure derivatives of distances, since the feature is 
already present in the distance, and this at smaller redshifts. Although this investigation 
certainly needs refinement, we already appreciate a signature of the different curvature 
evolution that furnishes a clearcut prediction for future observations (see [82] for details 
on the construction of template metrics, fits to observational data, and predictions). 

6.3. Strategies for generalizing numerical simulations 

The architecture of Newtonian simulations does not allow to describe the generic 
interaction between structure formation and the background geometry. This latter is, 
by the very nature of a Newtonian description within a coordinate frame comoving with 
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a (possibly relativistic) background solution, a foliation of non-dynamical Euclidean 
space sections. Moreover, kinematical backreaction terms vanish on the simulation box 
due to the technical restriction to a torus topology, see [IE] for the proof. 

Obviously, these restrictions have to be overcome, if we wish to conduct realizations 
of a physical cosmological model that go beyond the status of toy realizations to 
study structure formation. Thus, the need to construct relativistic simulations with 
dynamical geometry arises. Hereby one would wish that well-developed simulation 
techniques could be used and eventually interpreted within a relativistic framework. A 
straightforward method would be to integrate the ADM equations of general relativity. 
To date only special general relativistic systems are studied numerically, an effort 
to construct cosmological simulations awaits an attempt. In this situation we may 
ask whether one could improve the architecture of Newtonian simulations to render 
them relevant in a relativistic setting. A first example has been demonstrated in 
[125] , where the volume fraction between under-dense and over-dense regions has 
been used as approximate input into the relativistic framework. A measurement of 
fluctuations and the kinematical backreaction term could eventually be drawn from a 
Newtonian simulation and then, iteratively, taken as input into the average equations 
for a relativistic physical background. The introduction of comoving coordinates in such 
an evolving background will alter time-scales and the distance interpretation, but it is 
certainly a rough approximation and follows the spirit of using global template metrics 
as described above. 

A more systematic strategy is to first use analytical approximation schemes like 
the relativistic form of Zel'dovich's approximation, that is employed to initialize N- 
body codes. Work is in progress at this front and we shall discuss this latter strategy 
in concrete terms in the next section. 



7. Non— perturbative models for backreaction 

Backreaction models can be constructed on the basis of closure conditions on the 
averaged equations without specifying a metric (for discussions see [H] and [2T] j^. 
Examples include globally static and stationary cosmologies [21], scaling solutions 
[31], |125] . |113] . and specifications to other effective equations of state, e.g. |112] . 
These approaches only functionally depend on a metric and are to be considered as a 
motivated ansatz or as balance conditions imposed on the dynamics in the spirit of the 
virial theorem, where these conditions are not rooted back to explicit solutions of the 
inhomogeneous dynamics. It is here where a detailed investigation of inhomogeneities 
is needed to support or discard certain closure conditions. We already studied the 

I Note that the equations for averaged scalars form an infinite hierarchy of equations. Truncating them 
imphes the need for a closure condition. This happens in many physical systems, also for tensorial 
hierarchies, e.g. the velocity moment hierarchy of the Vlasov equation |30j . We can easily go down one 
level of this hierarchy, but always the need for closure conditions will arise, if the problem is restricted 
to a system of ordinary differential equations. 
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construction of homogeneous template metrics in the last section, which can be made 
compatible with global assumptions on the evolution of backreaction. Here we are going 
to study concrete inhomogeneous metric forms and the resulting backreaction models. 
We start with some general remarks. 

7.1. Some general notes on relativistic metrics and perturbation theories 

Consider for simplicity a spatial metric form g with coefficients Qij in an exact (co- 
tangential) basis dX^^dX^ . We can write any metric as a quadratic form of deformation 
one-forms, g = 6ab ® rj'', i.e. in terms of coefficients, gij = 6abV"'iV^^ Now, such a 
metric form is flat, i.e. its Ricci tensor vanishes everywhere, if there exist functions 
such that the one-forms can be written as exact forms, 77" = d/". In other words, if 
we can find a coordinate transformation = f°''"^{X\t) that transforms the Euclidean 
metric coefficients in a new basis, dx* ® dx^ , 5ijdx^dx^ = 5ahf°'\if^\jdX'^dX^, with a 
vertical slash denoting partial spatial derivatives, into the metric coefficients gij, then 
these latter are just a rewriting of the flat space. Given this remark, any perturbation 
theory that features metric forms of the integrable form, does not describe relativistic 
inhomogeneities; metric coefficients of the form g^j = Sahf°'\if^\j describe Newtonian 
(Lagrangian) perturbations on a flat background space. A truly relativistic perturbation 
theory deforms the background geometry; in other words, the perturbations live in a 
perturbed space, not on a reference background. This remark also shows that relativistic 
perturbation terms that are not related to coordinate artifacts can never contain full 
divergences, since this latter needs integrable one-form flelds. 

In light of these introductory remarks, an inhomogeneous relativistic metric 
produces curvature. The volume-average of this intrinsic curvature on some domain 
does not obey a conservation law, as can be explicitly seen in the coupling equation 
to the fluctuations (jS]). In particular, intrinsic curvature does in general not average 
out to zero; for details on curvature estimates see [29]. This fact in itself shows the 
existence of a dynamical evolution of an averaged curvature, as soon as structures form 
|18] . On the contrary, standard perturbation theory formulated on a flxed background is 
constructed such that the averages always vanish on the background, demonstrating the 
limited nature of results obtained by standard perturbation theory. Here, we identify the 
crucial difference between a fully relativistic cosmology and a quasi-Newtonian model: 
for the latter, the only fluctuating flelds are the matter flelds and for them we have a 
conservation law that assures that over- and under-densities compensate each other, 
even for a nonlinear density distribution, while for the former the curvature is also 
fluctuating; the argument that applies to the density distribution does not apply to the 
curvature distribution. 

Another perturbative argument aims to justify the validity of the homogeneous 
geometry, even down to the scales of neutron stars [72l [67] . As argued in [72] and 
by many others, perturbations of the metric remain small with respect to the flat 

§ We use indices i,j,k to denote coordinate indices, and the indices a,6, c as counters of e.g. forms. 
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background. While this is true, this does not contradict the existence of a large 
backreaction effect, since these latter depend on spatial derivatives of the metric that can 
be large |32], [78], |107] . Moreover, the perturbations are considered on a background 
that does not interact with structure. Perturbations may be small on a different 
(physical) background, in which case a perturbation may already live in a background 
with strong curvature (a zero-order effect). 

It is therefore not fruitful to argue against the relevance of backreaction within 
standard limited schemes, but rather an effort to generalize perturbation theory is 
needed. Efforts to construct a fluctuation theory around a physical average are the 
subject of current studies. Before such a more general theory can be constructed, it is 
necessary to first look at some results that hold for exact solutions and for approximation 
schemes describing truely inhomogeneous deformations as outlined above. We shall 
therefore discuss some results in the next subsection, obtained for the spherically 
symmetric LTB solution and for the relativistic Lagrangian perturbation theory. These 
results will be presented elsewhere within comprehensive investigations [Ml 135] . 

7.2. Backreaction for spherically symmetric solutions 

A large number of recent publications on the evaluation of backreaction is based on the 
spherically symmetric LTB solutions (for reference lists see the articles [581 13 EHl El I117[ 
[971 l93l 1119] . the book by Bolejko et al. [8], as well as the contributions [9] and [89] in 
this volume). A comprehensive study of LTB solutions in relation to the backreaction 
formalism discussed in this review and to the existing literature is provided by Sussman 
[119] . This paragraph just focuses on the special character of this class of solutions and 
provides some exact results. Before we come to the relativistic setting, we briefly recall 
the situation in the Newtonian theory [16], [25], |76j . 

Newton's Iron Spheres 

In Newtonian theory the background geometry is fixed and given by a Euclidean vector 
space, the Newtonian spacetime. Let the spatial domain of averaging V = be taken 
as a sphere with radius R. The velocity v inside is only depending on the distance r 
to the origin and always parallel to the radial unit vector e^, v = u(r) e^. By doing this 
we exclude rotational velocity fields. The chosen domain stays spherical at all times. 

The Newtonian velocity gradient is denoted by (wj ,,), with a comma indicating 
derivative with respect to a non-rotating Eulerian coordinate system; it may be 
characterized by its three principal scalar invariants, the trace I, the dispersion of its 
non-diagonal components II, and the determinant III (for explicit expressions we refer 
the reader to [52], [25]). The averaged first invariant may be obtained directly using 
Gauss' theorem: 

(I(^'<.)>., = 4^ X/'- V . V = ^ £^ dS . = 3^ , (17) 
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whereas the averages of the second and third invariants require some basic calculations. 
One obtains the relations |25] : 

= I , ; {mv^,))v, = ^ . (18) 

The first of the above relations implies that the backreaction term vanishes identically, 
Qsphcncai _ 2 (H)^^ — ^(^(Y)^^^ = 0, a result which is in accord with Newton's "Theorem 
of the Iron Spheres" . Since this relation holds true on every scale, the exact averaged 
equations (E]) reduce to the standard FLRW equations. 

Some exact results for LTB solutions 

The LTB solutions for dust generalize the well-known FLRW solutions for dust: their 
metric not only depends on the time-coordinate as in the FLRW model, but also on the 
radial coordinate. The spherical domain can be seen as a superposition of infinitesimally 
thick homogeneous shells governed by their own dynamics. In a comoving-synchronous 
setting (see e.g. [58] for a demonstration but with different notations) the line-element 
has the form: 

ds"" = -de + -5^!M^rfr2 + R\t,r)dn'' , (19) 
1 + 2E[r) 

E being a free intrinsic curvature function of r satisfying E{r) > —1/2; the prime 

denotes partial differentiation with respect to r. 

In this metric, the scalar parts of the Einstein field equations read: 

M being another free function of r related to the radial density profile; the overdot 
denotes partial time-derivative. Using the relation between the expansion tensor and 
the metric tensor in the coordinate form 6 ^ := ^g^'^gkj, the averaged scalar invariants 
of the expansion tensor on a simply-connected LTB-domain can be calculated |35] : 



4^ rrr, dr (i?/?^) 



4^ rr^ dr (R^R 



471- ri, dr [R^ 



^av^X '''' 



where the volume is given by 
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The deviation of the averaged scalar curvature from a constant-curvature model can 
also be averaged on a LTB domain: 

w,,,.-i5^r*ei*-6^25#. (25) 

VltbJo VTT2E Vl^l 

These integrals can be straightforwardly solved for E{r) = Eq = const, yielding the 
following exact results for the averaged scalar curvature, 

A{E{r)Ry AEq 

^^^^^^ = -^rT- = -WW) ' ^ ^ 

and for the averaged invariants: 

<"(0'.)>.TB = \ m))lTB ' (III(0^.)>.TB = ^ {m))lTB ' (27) 

Combining the averaged invariants into the backreaction term Qltb, cf. Eq. ([3]), we 
obtain for a spherically symmetric domain with a strong restriction on the curvature 
functioiljj]: 

Qltb = , Wltb = 0. (28) 

We here generalize to non-fiat domains a result obtained in [98]. Comparing Eq. fl25l) 
and the result for the averaged curvature, one can express k-^i as a function of Eq: 
kv, = -2Eo/R^{ti,rv). 

The result (!28l) mainly shows that spherically symmetric LTB solutions for a 
geometry with zero intrinsic curvature are quasi-Newtonian, i.e. they are too special 
and not useful to access the backreaction problem. Only work on LTB solutions that 
allow for a non-trivial curvature function E{r) and non-constant curvature geometries 
are relevant in this context. 

The result (|28l) can be interpreted as a generalization of what people have in mind 
when they quote Birkhoff 's theorem, since here the density distribution is continuous and 
Qltb = in general implies that the scale factor, volume-averaged over a spherically 
symmetric inhomogeneous distribution, follows the FLRW equations. Note, however, 
that it is by far not enough to have a spherically symmetric distribution, since the result 
(l28l) is very special and we cannot expect a similar theorem to hold in a more realistic 
situation. This point is important since, for this latter reason, we can expect to learn a 
lot from LTB solutions concerning the backreaction problem (see especially the review 
by Sussman [119] ). 

The above subclass of LTB solutions is contained in a wider class of backreaction 
models that are based on a relativistic Lagrangian approximation scheme. We are going 
to give some related results in the next subsection. 

II The restriction E = Eq = const, eorresponds to self-similar LTB solutions if we require at the same 
time that the function il/(r) oc r (R.A. Sussman, priv. comm.), see [117j . 
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7. 3. Backreaction in relativistic Lagrangian perturbation theory 

In this subsection we report on an application of relativistic Lagrangian perturbation 
theory for the construction of a generic backreaction model. The formulation of this 
theory as well as the formalism used to obtain the results sketched here will be published 
in forthcoming papers [MII35]. The idea of the construction of this theory is to employ 
as a single dynamical variable the spatial deformation one-form fields (Cartan's co- 
frames), 77", cast the full set of Einstein's equations for dust matter into a form that 
features only this variable, and then setup perturbation and solution schemes for this 
deformation field. Other fields like the 3— metric are then functionally expressed through 
solutions of this perturbation variable. 

We write the 3— metric as a quadratic form of non-normalized co-frames, 

gii = Gabr]\r,\ ; gi^it,) = Gi^{t-^ . (29) 

To choose non-normalized frames, that has been suggested by [IQ], bears the advantage 
that the resulting expressions are similar to their Newtonian counterparts (these latter 
can be found in [25]). 

The obtained perturbation solutions for the deformation one-forms read to leading 
order (and with a restriction of initial data that eliminates decaying modes): 

^^V.(i,^'):=a(t)(5'^, + e(t)^'^.) , (30) 

where ^^^^ = P'^iit-,, X^) , ^(t;) = 0, a{t{) = 1. The function ^(t) solves the well-known 
first-order equations to be found in [HI [151 [5]; RZA stands for "Relativistic Zel'dovich 
approximation", generalizing Zel'dovich's idea [131] and suggested first by Kasai [75] 
(for normalized co-frames). Contrary to Kasai's definition we consider the full 3— metric 
from first-order deformations: 

""^^.it, X'^) = a^t) {G,, + m {Gaj^% + G..^^) + e{t)Gab^%^'j} . (31) 
One then obtains for the RZA backreaction model in a non-normal basis: 

e (71 + + e'73) 

(l+^(I.),^+e2(IIi)c,+e(IIlOcj' ' 

where we have defined the set of initial data featuring the initial principal scalar 
invariants of the expansion tensor (the first is the initial backreaction term): 

72 :=6(IIIi)c, -^(IIi)c,(Ii)c, ; 

73:=2(Ii)c,(IIIi)c,-^(IIi)^^. (33) 

We note that for non-normalized co-frames the initial 3— metric tensor does not appear 
explicitly in the expression for the backreaction model; the domain V is, however, 
Lagrangian, i.e. it is frozen into the evolving metric. Initially, the domain V can here be 
chosen to be a section of a Euclidean space, denoted by Co- All relativistic expressions 
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have a straightforward Newtonian hmit by sending the deformation one-forms to exact 
forms 77" — )■ d/"=\ where the counting index for the forms a becomes a coordinate index 
i. Note that the Newtonian approximation contains integrable averaged invariants. 

In the situation of a spherically symmetric model restricted to the class of LTB 
solutions with trivial curvature function E{r) = Eo = const., the averaged invariants 
obey the (quasi-Newtonian) relations fl27|) . compatible with the backreaction model 
f p2|) . The generic model fl32|) therefore contains this class of LTB solutions in a subclass. 
Furthermore, the leading term in the model (l32l) agrees with the backreaction model of 
the linear perturbation mode, 



initial 



gUncar ^ ^2ginitial ^ ^ ^^^-^ 

(X 

where the latter equality holds for an Einstein-de Sitter background. This mode plays 
an important role in the evaluation of the backreaction effect, since it forms the weak- 
backreaction limit of an exact scaling solution [3T] : 

^initial yyinitial 
scaling = \ '^'^^""SW^, = , (35) 

that has to be compared with the competing sources in the balance equation of the 
averaged Hamiltonian constraint (the second of Eqs. ([2]))0, 

(^)p oc ^ ; kvo^\ ; ^'^^''"SQi, oc — ; A = const . (36) 

fl^j 0"£) Old 

In Section |5] we discussed a model that assumes the scaling laws ( l35l) on subdomains 
of a multiscale cosmology where the global evolution mimics a cosmological constant 
behavior as a result of the expansion variance between the subdomains. 



l.Ji^. Backreaction models for relativistic scalar metric inhomogeneities 

Thus far we worked with the matter model 'irrotational dust'. To discuss more general 
matter models such as radiation or the fluid picture of a scalar field, we have to briefly 
recall a covariantly defined set of averaged ADM equations obtained previously |19j . 
These equations are valid for any spacetime foliation within the class of foliations with 
vanishing shift, and for any choice of the lapse function and the inhomogeneous 3- 
metric. We shall then specify the 3-metric, investigate scalar metric inhomogeneities 
for any choice of the lapse function, and evaluate the relevant backreaction terms. This 
we can do in general without resorting to any other approximations than those implied 
by the restriction to scalars and irrotational flows. We shall also discuss the evolution of 
backreaction in the longitudinal gauge, and we shall put forward crucial arguments in 
favor of a non-perturbative versus a perturbative interpretation of backreaction using 
this example. 

% For some further reniarks in this context see [53], Sect. 4.2. 
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7.4- 1- The averaged ADM equations for vanishing shift In this subsection we entirely 
follow the notations and results given in [19]. We shall study spatial averages in a 
hypersurface defined by the choice of the in general inhomogeneous lapse function 
N{X^,t) and inhomogeneous 3— metric coefficients gij{X'^,t) in the line-element 

ds^ = -N^dt^ + gijdX'dX^ . (37) 

(X* are local coordinates in a t = const, hypersurface). This line-element is sufficient 
to analyze the example of the so-called longitudinal gauge that we shall consider below. 

Perfect fluid sources are characterized by a diagonal energy-momentum tensor with 
energy density e and pressure p, T^y = eu^u^+ph^^. We may choose to project onto the 
fluid's restframe, deflning the projection tensor through /i^jj, = g^^ + u^Uy, i.e. we project 
onto hypersurfaces orthogonal to the fluid's 4-velocity u^. We employ the 4-velocity of 
the flow in the form 

; ft = £±Z , (38) 

n g 

together with the decomposition into kinematical parts of the 4-velocity gradient, 

u^^-u = \oh^y + a^y + Uf,y - ix^Uy , (39) 

where the inhomogeneous normalization of the 4-velocity gradient h is given by the 
injection energy per fluid element and unit restmass, de = hdg with the restmass density 
Q [13]; is the rate of expansion, a^y and uj^y the shear and vorticity tensors, respectively. 

The existence of a scalar 4-velocity potential S together with the choice f l38l) implies 
that the conservation equations T^^'^.^ = are satisfled, but also that the flow has to be 
irrotational and that the covariant spatial gradient of S (denoted by a double vertical 
slash in this paper) vanishes [121 [131 EO], [19] : 

= h^"hfu[a;p] = -h^^hf ( d[^S/h ) .^j = ; = h^^d^S = d^S + u^S = , (40) 

with the covariant time-derivative S := u^S-^ = h. For the special case of an equation 
of state of the form p = 75 we obtain: 

e=^h'^'^^ ■ p = h'+'/'' ■ f, = l±2/,i/7 . h + ^9h = ; h ^ S , (41) 

which reduce to the familiar expressions for a free minimally coupled scalar field source 
(a "stiff fluid" with 7 = 1; this case has been exploited in [26]). 

The averaging operation in terms of Riemannian volume integration is performed, 
as in the dust case, over the hypersurfaces orthogonal to u^, restricting again attention 
to the scalar functions \&(X*,t), cf. Eq. ([1]). We also consider the same deflnition as for 
the dust case of a dimensionless volume scale factor aT>(t), which implies that we are 
only interested in the volume dynamics of the domain; is a functional of the domain's 
shape (dictated by the metric) and position. As in the dust case we require the domains 
to follow the flow lines, so that the total restmass Mp := jj^gJd^X contained in a 
given domain is conserved. With a non-constant lapse function we have, however, to 
introduce a scaled (t-)expansion 6 := NO, which describes the rate of change of the 
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domain's volume expansion in the spatial hypersurfaces, that on average defines an 
effective Hubble-functiona0: 

= ^ = 3*^H . 3^^^ =: 3i?. . (42) 

Vuytj ax> ax> 

For an arbitrary scalar field T(X*,t) we make use of the non-commutativity relation: 

(T)^-(T%= (T^)2,-(T)^(^)i,, or, alternatively, (T)^ + 3/7p(T)p = (r + T^)^.(43) 

Averaging Raychaudhuri's equation and the Hamiltonian constraint, we can cast 
the resulting equations into a compact form (to be found in Corollary 2 in |19j): 

3^ + 47rG (£efr + 3pefr) = ; 6//|, - IQ-kGe^^ = ; e',s + 3Hj, {e,s + Pcs) = , (44) 

with the following fiuctuating sources: 

levrG^eff := lQ7iG{e)v - Qv - {n)v , (45) 

167rGpeff := 167rG(p)2, - Qv + - ^Vv ; (46) 

e := N'^e and p := N^p are the scaled energy density and pressure of matter, respectively. 
The kinematical backreaction term is given by: 

Qr,:=2{N'lI)v-l{Ne)l ; (47) 

it is built from the principal scalar invariants 211 := 6"^ — K^jK\ and K^^ = —9 of the 
extrinsic curvature, K'^^ = ~\9^^j^g'k'j- The averaged 3— Ricci scalar curvature TZ and 
the acceleration terms {dynamical backreaction) read: 

{n)v ■■= {N^n)v ; Vv ■■= {A)v + {^^)^ . (48) 
with the scaled (t—) acceleration divergence A := N'^A = A^A^'|||f]. 

7.4-2. Scalar metric inhomogeneities in a metric form corresponding to the conformal 
Newtonian gauge The so-called conformal Newtonian or longitudinal gauge is often 
employed in the study of perturbations on a Friedmannian background cosmology, and 
is considered a preferred frame because it offers a well-defined Newtonian limit |124j^l . 
The topic addressed in this review, i.e. the impact of inhomogeneities on expansion 
properties of the Universe (backreaction) is also often discussed in this gauge [61], [62] . 
[96| [1], [651 ESI [90], [87], and many others. In the recent paper [59] it is shown, in an 
infiationary model, how the observers in this gauge are related to the free-falling ones. 
Even a 'no-go conjecture' has been raised on the issue of whether backreaction can be 

+ Wc shall reserve the overdot for the covariant time-derivative (defined through the 4- velocity u^): 

:= = jf§[, and we abbreviate the coordinate time-derivative by a prime in the sequel. 

* A single slash denotes partial differentiation with respect to the coordinates X*, and a double vertical 
slash covariant spatial differentiation with respect to the 3-metric as before. 

H Note that the framework discussed in Subsection 7.3 also offers a well-defined Newtonian limit. 
However, there one falls on the Lagrangian form of the Newtonian equations |51[ 153) . 
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significant [72], also advocating the post-Newtonian metric as a sound model for most 
cosmological studies. George Ellis gives a related discussion in this volume [55] . 

Since the averaged equations briefly reviewed above are valid for any choice of the 
lapse function and any ansatz for the 3— metric, we are in the position to calculate 
backreaction effects with some generality, i.e. with no need to invoke approximations 
other than those implied by the restriction to scalars and irrotational flows. We shall 
do this calculation explicitly with the aim to illustrate the strict, non-perturbative 
application of the post-Newtonian metric form as a solution of general relativity, and 
to learn some issues about common perturbative interpretations of this metric form. 

In the longitudinal gauge the lapse function and the 3— metric are specified in many 
studies as to provide a "Newtonianly perturbed model" in the following form: 

iV2 = i + 20 ; = (1-2^) 7,, , (49) 

with the scale-factor a of a homogeneous-isotropic background model, and a constant- 
curvature 3— metric 7jj. For simplicity, we are going to choose the Euclidean metric 
lij ~ ^ij what follows. From what has been said previously we are in the position 
to evaluate all the variables, in particular the backreaction terms, as functionals of the 
lapse function N{X\ t), the metrical inhomogeneities ip{X\ t), and their time and space 
derivatives. 

Note already the subtle element that in general relativity X* are to be local 
coordinates in a perturbed space. In turn, factoring out a scale factor as in the metric 
form (H9|l implies that, only if = '0 = 0, the scale factor obeys the standard Friedmann 
equations with respect to the coordinate-time; in the perturbed space the scale factor 
acquires a dependency on X* as seen in the hypersurfaces t = const., defined by the 
inhomogeneous lapse function. In the presence of perturbations it makes no physical 
sense to factor out a function of t that obeys Friedmann's equations. If the background is 
not perturbed, then the lapse function can only be time-dependent (for the background 
equations in hypersurfaces with a time-dependent lapse see [19], Sect. 4.1). 

Note also that the simple ansatz for the (conformally flat, i.e. vanishing Cotton- 
York tensor) metric requires that a scalar function models all six metric components. 
It is therefore expected that in the case of a strict application of this metric form as 
a solution to general relativity we are dealing with a highly restricted situation. We 
remark that this metric ansatz is general for the case where the trace-free symmetric 
part of the extrinsic curvature (the shear tensor aij) vanishes: 



1 N 



^' a' 



(5V ^ Ry - Ia'V^- =: = . (50) 

However, the post-Newtonian form is employed with the implicit understanding that 
101 and are small compared to 1 (together with corresponding requirements on their 
derivatives |124] ). For perfect fluids and small peculiar- velocities (j) = tp. The so- 
constructed model is designed to be in a "near-Friedmannian state". In what follows 
we have to keep in mind that Eq. (!50|) implies that any application of this metric-form to 
describe inhomogeneities can only be considered in an approximate sense. Especially for 
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a non-tilted slicing where the expansion tensor is proportional to the extrinsic curvature 
(peculiar velocities are not small but vanish), we have vanishing shear everywhere and 
this implies homogeneity in cosmologically relevant case iffl . 

In the sequel we shall consider ( 149|) as an ansatz for the metric and only later, in the 
final result, we may look at the importance of the various terms; we do not invoke any 
further approximation in the following calculations. The usual practice of employing 
perturbative assumptions from the outset may mask the simplicity of the problem. Also, 
we prefer to not specify the lapse function; we shall retain so that we can discuss 
the result for different foliations. We introduce, to shorten the notation, the auxiliary 
variable with its derivatives with respect to the coordinate time: 

a:^-lnV^ ; ^ ; ^ ^ + . (51) 

Kinematics of the volume Let us start with the simple observation that, given the 
3— metric in the form (l49l) . the kinematics of the volume is determined by a given 
solution for ip: the volume of an averaging domain, and hence the effective scale-factor 
a©, is calculated from its definition Vt> = JcPX with J = ^J(^.ei{gij) = a^(l — 2^/^)'^/^: 



Vv = a^Vvi= / a''[l-2'ipY'' d-'X . (52) 
Jv 

The rate of volume expansion of a domain in the spatial hypersurface, written in terms 
of the Hubble functional Hx>(t), reads: 

\{Ne)v = \^ = H^= {H)v - («% , (53) 

with the local Hubble function H{X\t) := a' /a of the background model as seen 
in the spatial hypersurfaces specified by the inhomogeneous lapse, i.e. H = a/ a = 
a'/{aN) = H/N. At this place note that there is no ambiguity concerning the notion 
of averaged volume expansion, once the lapse function has specified the foliation of 
spacetime. Although being unambiguous, we have to come back to this point later, 
since a majority of papers on averaging scalar metric inhomogeneities employs Euclidean 
volume averaging on the background metric in a frame of global coordinates, which do 

ft Note that a priori the tilt of the 4— velocity relative to the hypersurface normal is not specified. For 
vanishing tilt, as considered here, and in the case of dust matter, shear-free motion implies homogeneity; 
this also holds true for large classes of perfect fluid models, see |46[ l47l |45] . We briefiy show for the 
case of dust matter that we can determine the lapse function such that the model is hypersurface- 
homogeneous: we use the momentum constraints, K^j^^ - — K'^^j = 0, for the extrinsic curvature of 
Eq. ((50)) . and integrate them to yield Ns{t) = ip' /{I — 2-0) — H, with a time-dependent function of 
integration s{t) that reflects the freedom of time-reparamctrization. On the other hand wc have from 
S = N~^S,t = h, with h ~ (e +p)/g: N — S,t/h. Equating the two relations for the lapse function 
gives ip' /{I — 2ip) — a' /a — s{t)S{t)/h. This shows that, if we require p = and hence N^i = and 
h = 1, the function can only be time-dependent and does not describe perturbations in the considered 
hypersurfaces; e.g. the Ricci tensor, given further below in explicit form, vanishes everywhere, since the 
space gradients of ip{t) vanish. For further discussions of this metric form including estimates, see [32| . 
Work is in progress that analyzes the present issue in a tilted slicing and non- vanishing shift vector. 
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not exist in a general relativistic setting, if the space is deformed by inhomogeneous 
perturbations (see the related discussion of George EUis [55]). 

The second time-derivative provides the kinematical equation for the rate of volume 
acceleration: writing 



a" 



-^ = H^ + Hl, (54) 

and, using the commutation rule (H3|) for both terms on the r.-h.-s., we calculate: = 

{H% - {a'% + Sliia'Yh - («%] - 6[{Ha% - {H^ia'U + 3[{H')^ - {H)l], and 

H-^ = {a')'^ — 2{H)D{a')T) + {H)v- Summing up the above terms, using H' + H'^ = a"/a, 
and rearranging some terms, we get for the rate of volume acceleration: 

^ = (-)+2[ ((«')%-(«%] 
a-D \ a / V 

-4[ {Ha')v - {H)v{a')v] + 2[ {H')t, - {H)l] 

~{a" + 2Ha -{a'f)v , (55) 

where a" /a = N'^{a/a) + N\a/a) is the local rate of acceleration of the background 
model as seen in the spatial hyp ersurf aces. 

This latter equation is quite simple given the fact that the second line above only 
appears, since we have factored out a "background model" which, within the spatial 
hypersurfaces, appears as a fluctuating background: there is a fluctuating t— Hubble 
function, which with the replacement H = NH results in a frame backreaction term. 
This term has to be taken seriously, since there are no "background observers" from 
which we could see a homogeneous behavior of the background expansion. 

Given a solution of Einstein's equations for ip and given a lapse, e.g., as a functional 
of '0, the above equation provides the answer to the question, whether the universe model 
'accelerates'. We did not use any dynamical equations of general relativity so far. 

Backreaction terms Since we did not employ any approximation beyond those of the 
chosen framework, the first and second t— derivatives of the scale factor calculated above 
are in accord with those found from the general equations (14^ . To demonstrate this we 
have to employ dynamical equations of general relativity. We evaluate the kinematical 
and dynamical backreaction terms: 

% = [ ((«')')^ - - 2[ {Ha% - {HMa')v] + [ {H')v - {H)l] ; (56) 

'^v = {A)r> + 3{^{H-a'))^ . (57) 

We employ the Hamiltonian constraint and Raychaudhuri's equation in the forms: 

n = IQnGe + 2AN^ - QH'^ + 12Ha' - 6(a')^ (58) 

a" + 2Ha' - {a'f = ^{e + ?>p)- ^{H -a')~l-A+-. (59) 

6 JM 6 a 
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Then, we take spatial averages, use the averaged t— scalar curvature {N'^TZ)d, and insert 
the average of the left-hand-side of fl59|) into the last line of Eq. fl55l) : we so find 
consistency with the averaged equations (144|) : 



— = -^{e + mv + UQv + Vv) ; GH^ = 16nG{e)r, - Qv - {n)v ■ (60) 
(The t— averaged curvature term cancels in the averaged Raychaudhuri equation.) 

7.4-3. Dynamical equations for scalar metric perturbations The Hamiltonian constraint 
and Raychaudhuri's equation may be rewritten by replacing the coordinate time- 
derivative (denoted by a prime in the previous equations) through the covariant 
derivative (denoted by an overdot), and use of the relation H = HN: 

n = IGnGe + 2A - 6H^ + 12Ha - 6a^ ; (61) 

a + 2Ha-a^ -^^{5e + 35p) + ^A = , (62) 

where we have split off the background model through introduction of the deviations 
5e := e — Eh and 5p := p — ph, and where the background quantities have to obey the 
covariant equation 

a AttG , , a!' AnG , , N' ~ 

- + ^{eH + 3pH) = ^ — + ^{£h + 3ph) = ^H. (63) 

Qj tj Qi O iV 

Equation fl62|) is key to the determination of solutions for tp; we therefore write it in 
terms of ip: 

- • 1 tij^ 

il: + 2Hi:- -{4nG[Se + 35p]~A){l-2ij)+ ^ =0 . (64) 
3 1 

The acceleration divergence A = A/N"^ in these equations can be explicitly expressed 
through ip only, if the lapse function is specified, e.g. through ip. If it is specified 
according to fH91) we obtain: 

IV 1 + 20 (1 + 20)^ 

involving the Laplace-Beltrami operator on the 3-metric, (note: F*.^ = (In J)|fc and 
J := v/d^t(^ = 03(1-2^^)3/2): 

A,0 := cp\ = = [^?^Vb-]|. + /Vb-r^fc = JW'J^\A\^ ■ (66) 

We notice that we have to specify the form of the lapse function. Then, we can proceed 
by solving or approximating the local evolution equation governing ip. Also this problem, 
if approached perturbatively from the outset, could mask the transparency of the general 
problem, and so may give rise to ambiguities in interpretation. 



-4=+^ = T^-7T^ . (65) 
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Example As an example that connects the above thoughts with the model in 
Subection 17.31 we shall specify the matter model to an irrotational dust continuum, 
i.e. p = 0, and therefore /i = 1, and e = q. We consider the comoving synchronous 
gauge with = 1 without loss of generality. Hence, we have no acceleration in the 
hypersurfaces, ^ = 0, and retain the simplified equation: 

^ + 2H^- ^6gil - 2^) + = ; 5g = g-gH=:gHS , (67) 

3 1 — Zip 

where the overdot is, in this gauge, equivalent to the coordinate time-derivative, and S 

is defined as the conventional density contrast. Since the energy density is reduced to 

the restmass density g, we can employ the general integral 

T_l . 1 Qinitial 1 ~l~ ^initial ^cq\ 

Q-Q^.ua,J ^.e. ^'-^(i_2^)3/2 2^)3/2 • (68) 



Although a general solution to Eq. f l67jl may exist and could be found, we here resort to 
approximations. Linearizing the above equation with respect to ip, i.e. also expanding 
5^(1 - 2iIj) to linear order, Sg{l - 2^)) ^ ^^[(1 + Sinitiai){l + Sip) - 1](1 - 2iP) ^ 
gniSinitiai + 3'(/'), we obtain a familiar equation for the approximate solution ip^ [85] : 

+ 2H1P'' - AnGgn^P'' = . (69) 

3 

Interestingly, this approximation may be interpreted as the scalar part of the (first- 
order) relativistic Lagrangian approximation discussed in Subsection 17.31 For, if we 
write the 3-metric in terms of Cartan one-forms, Qij = SabVlV^'j = a^(l — 2ip)6ij, the 
deformation is given by rjl = ay/1 — 2ip 5". On the other hand, the equation for the 
trace of first-order Lagrangian perturbations reads |3l] : 

^n-. = a[5\ + P'';\ ■ P:=P\ = 5akP\ ; P + 2HP - AT^GguP = A^TGgH5,nit^al. (70) 
Linearizing the trace of the Cartan deformation above we obtain the relation P = Sip"^, 
so that Equation ( !70|) implies Equation fl69|) . Keeping only the growing part of the 
homogeneous solution we have: 

= i^^nu^al{X')m with ^ + 2Hi - AttG gni^ + 1) = , (71) 

with the well-known solutions C,{t) of the linear perturbation theory, or the standard 
Zel'dovich approximation [13H [T^ [T5| [5] . 

Inserting this approximate solution into the general averaged equations results 
in a non~perturbative backreaction model. The averaged intrinsic curvature can be 
evaluated in two ways: firstly, from averaging the Hamiltonian constraint (l58l) as above 
(which is enough for the kinematically averaged equations), and, secondly, from the 
geometry of the spatial hypersurface: locally, we have for the 3— Ricci tensor: 



TV ■ = — 



(1 - 2i\)Y (1 - 2il)) 

and for the scalar 3— curvature Tl = TZ'' 



(72) 



a'' 



M , 3 {ViPf ^^3^ 



_{l-2ipy 2(l-2^)3_ 
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with the Euchdean operators Aip = S'^^ip^n^, Vip = ip\k- (Note that, since = 1, the 
scale-factor a depends only on the coordinate-time t; for non-constant lapse the above 
expression also contains spatial derivatives of this scale-factor.) Averaging the above 
curvature expression can be used to test consistency and therefore the quality of the 
proposed non-perturbative approximation for the backreaction model. 

Discussion: non-perturbative versus perturbative approaches 

The above example calculation for metrical inhomogeneities of dust furnishes a 
perturbative approach with regard to the local solution. Considering the backreaction 
problem, i.e. the evolution of average properties given the local solution, we are entitled 
to look at the averaged equations in their general form and only approximate the local 
evolution of ip, and this as general as we can (compare here some results obtained in the 
Newtonian theory [22]). In such an approach no approximation is done on functionals of 
ip like the averaged equations, the averaging operator itself, etc. We so deal with a non- 
perturbative result. The perturbative result should arise as a limit; however, this limit 
process is drastic: it involves not only linearizing the equations, but also the averaging 
operator, which explains that in previous work many more backreaction terms were 
found [96] ; compare also [12], I122j and the explicit discussion of this point in [13] and the 
contribution [S] to this volume. Exploiting our result in this limit is of course possible, 
but from a physical point of view not necessary. The present suggestion of a non- 
perturbative approach consists in performing approximations only for the fluctuating 
local sources, but how these approximations enter functional expressions like the volume 
average and the backreaction model is not further approximated. 

The discussion of scalar metric inhomogeneities also should make clear that for 
a realistic modeling of inhomogeneities the tensorial degrees of freedom have to be 
taken into account. As our example and its connection to the relativistic Zel'dovich 
approximation show, this can be achieved and forms the subject of forthcoming 
investigations [Ml EH] . 
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